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Abstract 

We investigate numerically the spin polarization of the current in the presence of Rashba spin- 
orbit interaction in a T-shaped conductor proposed by A. A. Kiselev and K.W. Kim (Appl. Phys. 
Lett. 78 775 (2001)). The recursive Green function method is used to calculate the three terminal 
spin dependent transmission probabilities. We focus on single-channel transport and show that 
the spin polarization becomes nearly 100 % with a conductance close to e 2 jh for sufficiently strong 
spin-orbit coupling. This is interpreted by the fact that electrons with opposite spin states are 
deflected into an opposite terminal by the spin dependent Lorentz force. The influence of the 
disorder on the predicted effect is also discussed. Cases for multi-channel transport are studied in 
connection with experiments. 

PACS numbers: 
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I. INTRODUCTION 



Recently, spin-dependent electronic transport is attracting considerable attention because 
of possible applications to spintronics Q, Q|. Many of the proposals for two dimensional 
(2D) spintronic devices are based on the presence of spin-orbit coupling in the 2D electron 
system (2DES) semiconductor heterostructure. There are two types of spin-orbit coupling 
terms in such systems. One is the so-called Dresselhaus term which originates from the 
inversion asymmetry of the zinc-blende structure [3| . The other is described by the Rashba 
Hamiltonian, 

OL 

71-R = -r{°xPy - °yPx) (1) 

where a denotes the strength of spin-orbit coupling, a, and Pi (i = x, y) are the Pauli 
matrices and the components of the momentum, respectively. The Rashba mechanism is 
due to the effective electric field originating from the asymmetry of the potential confining 
the 2DES 

It is well known that the Rashba term dominates in narrow-gap semiconductors while 
the Dresselhaus term is dominant in wide-gap systems jf]]. Since the strength of the Rashba 
term can be controlled via spin-orbit interaction 

have become most promising for spintronic applications. 

In order to realize such devices, one needs spin polarized electrons in the semiconductor 
inversion layer. Most straightforwardly, one could generate spin polarized electrons by at- 
taching ferromagnetic metallic contacts to the 2DES and by injecting a current However, 
it has been found that in practice the efficiency of the spin injection from a ferromagnet into 
a semiconductor is very poor because of the conductivity mismatch [10]. Thus, alternative 
methods have to be invented. Since strong spin-orbit scattering can lead to spatially in- 
homogeneous spin polarization ^l| the generation of spin polarized electrons via spin-orbit 
coupling is in principle possible. 

In this paper we consider a conductor with a T- shape structure with Rashba spin-orbit 
coupling as originally proposed by Kiselev and Kim .12]. For relatively small strength of 
the spin-orbit coupling, these authors have obtained high spin polarization. However, the 
corresponding conductance has been found to be very small. This problem of the small 
conductance has been eventually overcome by considering a ring-shape electron resonator 
j"!^ . It has also been reported that spin accumulation occurs for considerably strong spin- 
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orbit coupling in a quasi-lD wire 

In both of the above cases, electron transport in the lowest subband (single-channel 
transport) has mainly been considered. Due to the self-duality of scattering matrix j3| for 
the system with spin-orbit interaction, no spin polarization of the current can be obtained 
for single-channel transport in two-terminal devices. Therefore, one has to consider at least 
a conductor with three terminals. The single-channel limit seems to be ideal because one 
can completely suppress the effect of the D ' yakonov- PereP relaxation, which is the major 
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spin relaxation mechanism in such systems [16 , 

We will show in this paper that the amplitude of the spin polarization becomes almost 
perfect with only very little loss of the conductance if the spin-orbit interaction is sufficiently 
strong. We argue that the predicted effect should be experimentally accessible in InAs. 

In the next section, we describe the model system to be investigated numerically. We 
use the Ando tight-binding Hamiltonian with Rashba coupling in the off-diagonal matrix 
elements j^J. In Sec. IIH[ the results for the dependence of the polarization on the energy 
and the strength of spin-orbit coupling are presented. The amplitude of spin polarization 
is shown to depend on the ratio between the 7r-phase spin precession length and the width 
of the quantum wires of the T-shape conductor. In Sec. IIVI we discuss the origin of 
the polarization by investigating the spin states of the wave function. We show that the 
propagating electrons are deflected at the junction by "Lorentz force" due to the spin-orbit 
induced effective magnetic field proportional to the z-component of the spin state. The 
effects of disorder and other channels are also investigated. The final section is devoted to 
the summary of this paper and the discussion for experimental realization. 



II. THE MODEL 

We consider the T-shape conductor shown in Fig^ in the presence of Rashba spin-orbit 
coupling. The sample region is connected to three electron reservoirs by ideal leads. Electron 
current is injected into the sample from the reservoir 1 and goes to reservoirs 2 or 3. At 
small voltages, the currents I21 and J31 from reservoir 1 to reservoirs 2 and 3, respectively, 
are proportional to the conductances G 21 and G 31 . 

The Rashba spin-orbit coupling can be described in the tight-binding language by the 
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FIG. 1: Schematic view of the T-shaped conductor. Current injected from reservoir 1 can go to 
reservoirs 2 or 3. Shaded: regions with non-vanishing spin-orbit coupling; parameters: N w = 10a 
and Ni = 20a with a lattice spacing of tight-binding model. 



Ando Hamiltonian 



where 



H = Y^ Wic\ a c i(7 - v i°,jv' c lcjv> , ( 2 ) 

i,cr {i,'j),cr,a' 



, cos 9 — sin 9 . 
Vi, i+& = V | | , (3) 

sin 9 cos 9 



and 



, cos 9 i sin 9 , 

Vi,i+y = V | | . (4) 

i sin 9 cos 9 

Here, c\ a (ci a ) denotes the creation (annihilation) operator of an electron on site i with 
spin o", Wi the random potential on the site i distributed uniformly in [— W/2,W/2], and 
Vi,i+x(Vi,i+y) the hopping matrix elements in x-{y-) directions. The hopping is restricted 
to nearest neighbours. The hopping energy Vq = h 2 /2m*a 2 , where m* is the effective 
electron mass and a the tight-binding lattice spacing, is taken as the unit of the energy. 
The parameter 9 represents the strength of the spin-orbit coupling and is related to Rashba 
coupling a by 

a~29V a (for0<l). (5) 

The conductance between reservoirs J to I J^jJ and the corresponding spin polarization 
are defined by 

G 7J = G Tr4 J t 7J , (6) 
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and 

= ~ 4- (k = x,y,z), (7) 

with the conductance quantum Go = e 2 //i, t/j denoting the transmission matrix from reser- 
voirs J to /. Below, we will focus on the transport between reservoirs 1 and 2. Transport 
between reservoirs 1 and 3 can trivially be deduced via current conservation and the sym- 
metry of the system. We calculate the amplitude of the total spin polarization defined by 

\P\ = {P 2 X +P 2 y +P 2 z ) l '\ (8) 

instead of considering only the ^-component. 

As described in detail in Appendix El we can obtain the transmission coefficient for 
the incident channel v with velocity v „ in the probe J and out-going channel fi with velocity 
Vfj, in the probe I as j^] 

C = (^) [—(U I )~ 1 G IJ U J {(A J )~ 1 - A J }U (9) 



where G IJ is Green function corresponding to the transport from the probe J to I. A 
contains only diagonal elements, A J (i,j) = Xf ■ Sy, where A/ is the z-th eigenvalue of 
the transfer matrix for an ideal region in the probe J and U 1 ^ consists of the set of 
eigenfunctions for {A[ ( ' J ^}. 



III. RESULTS 



Using the recursive Green function method, we calculate the amplitude of the spin po- 
larization. For the system size we assume N w = 10a and iVj = 20a. We first consider an 
clean system (W = 0). 

Figure 121 shows the dependence of the conductance and the spin polarization on the Fermi 
energy E for weak and strong spin-orbit couplings. We consider the energy region for single- 
channel transport. For weaker spin-orbit coupling (6 = 0.027r), high spin polarization is 
obtained for energies just before the second channel opens (E ~ — 3.68Vo). The correspond- 
ing conductance is small as compared to Go [12]. Almost perfect polarization is obtained 
together with a conductance close to G for stronger spin-orbit coupling {6 = 0.067r). Here, 
the polarization is almost insensitive to the energy except near the band edge. 
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FIG. 2: Conductance G and spin polarization \P\ as functions of the Fermi energy E in the 
region of single-channel transport. Almost 100% polarization is obtained together with G ~ Go 
for stronger spin-orbit coupling (9 = 0.067r). For 9 = 0.02-7T high polarization is obtained only at 
energies just before the second channel opens. 

Figure |3] shows the dependence of the conductance and the spin polarization on the 
strength of spin-orbit coupling at energy E = — 3.8Vo. With the increase of the strength of 
spin-orbit coupling, P y also increases monotonicaliy while P x and P z oscillate. We also note 
that the conductance increases together with the amplitude of the polarization. Due to the 
symmetry of the T-shaped conductor, the current into reservoir 3 has the same polarization 



in the direction of y but opposite polarizations in x- and z-directions 
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FIG. 3: Conductance G and spin polarization P^ in the directions of A: = x,y,z as functions of 
the strength of spin-orbit coupling 9 for E = — 3.8Vo; P y increases monotonicaliy with increasing 9 
while P r and P z oscillate. 



What are the conditions for achieving almost perfect polarization? We define the 7i~-phase 
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spin precession length L BO (\P\, N w ) = ira/29(\P\, N w ) where 6*(|P|, N w ) is the spin-orbit cou- 
pling strength giving rise to the polarization \P\ for the width N w . From the plot L so (|P|, N w ) 
as a function of the width of the system N w (Fig. HJ) one concludes that L so (\P\, N w ) is al- 
most linear in N w and high polarization, \P\ > 0.99, is achieved for L so (\P\, N w ) < N w . On 
the other hand, no dependence on the length of the leads (Ni) is observed as shown in FigO 




5 10 15 20 
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FIG. 4: The spin precession lengths L so for \P\ = 0.5,0.75,0.99 as functions of the width of the 
wires N w . Energy is set to be the middle of lowest and first excited subband, e.g., E = — 3.8Vo 
for N w = 10. The polarization becomes almost perfect if the spin precession length becomes 
shorter than the wire width, L S0 (|P|,A^) = 7ra/29(\P\, N w ) < N w . Solid line: L so (\P\,N w ) = N w 
corresponding to \P\ ~ 0.97. 
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FIG. 5: Conductance G and spin polarization |P| as functions of the wire length ./V; for several 
strengths of the spin-orbit coupling at E = — 3.8Vo- The polarizations are insensitive to the length 
of leads except for very short length. 
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IV. DISCUSSION 



It is well known that the spin-orbit interaction is the origin of the anomalous Hall effect 
in ferromagnetic materials [23] . This is because the spin-orbit coupling affects polarized con- 
duction electrons as a spin- dependent effective magnetic field. In this section, we investigate 
the spin wave function and show that propagating electrons are deflected at the junction 



by the Lorentz force due to the 



the ^-component of the spin 
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spin-orbit induced effective magnetic field proportional to 



24. 



251 ] . We also investigate the effect of disorder and the 



influence of other transport channels. 



A. Spin-orbit induced effective magnetic field 

A convenient way of deriving the effective magnetic field is to estimate the flux per 
plaquette from the Aharonov-Bohm phase. Let an electron initially be at When it 

hops from to (i + x,j) to (i + x, j + y) according to Eqs. J3J) and (HJ), the SU(2) phase 
exp(— id<j x ) exp(i8a y ) is acquired. On the other hand, when it hops from (i, j) to (i, j + y) 
to {i + x,j + y), the phase it obtains is exp(i9a y ) exp(— \9a x ). Therefore the interference 
between the two path along the plaquette becomes 

(exp(i9a y ) exp(— i#<7 x ))* exp(— \9a x ) exp(i9a y ) = exp(i9a x ) exp(—i9a y ) exp(— \9a x ) exp(i9a y ) . 

(10) 

This can be evaluated on the basis of Campbell-Hausdorff formula, 

exp(^X) exp(flr) = exp (e(X + Y) + y [X, Y) + 0(9 3 )^j , (f f ) 

to find 

exp(i9a x ) exp(—i9a y ) exp(— \9a x ) exp(i6'cr y ) (12) 
w exp(i(a x . - ay) + \9 2 a z ) exp(-i9cr x ) exp(i9a y ) 
w exp(i(cr 1 . -<J y - a x ) + 2i9 2 a z ) exp(i9a y ) 
w exp(2i^V 2 ) , 

where we have dropped terms higher than 9 2 . Using the relation, 

27iBa 2 /(h/e) = 29 2 a z , (13) 



we have 124, 2 
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H29 2 



e a 



2 °* 



(14) 



We now discuss the condition for the high polarization in terms of this spin dependent 
Lorentz force. We use Eq. (jl4j) to estimate the strength of magnetic field 



m - A - h - ( B - 

7r e V a 



(15) 



where we have averaged o z to be 2/7T, since the variation of the expectation value of a z is 
described by the cos-function. Since the kinetic energy is comparable to the confinement 
energy in single-channel transport, the velocity of an injected electron can be assumed to 
be 

(16) 



H 7T 

m* N w 

Then the corresponding cyclotron diameter is given by 



2l B {9) 



2m* v 212 



(17) 



eB{6) N w 

The cyclotron diameter becomes shorter than the wire width (21 §(6) < N w ) if the spin 
precession length becomes shorter than the wire width (L so < N w / v2). As a result, electrons 
with opposite ^-component spin are almost completely separated at the junction and nearly 
perfect spin polarization is obtained (FigEJ). This situation is similar to the mesoscopic cross 
junction in magnetic fields 
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FIG. 6: Schematic view of the electron trajectory. A pair of the electrons with anti-parallel spins 
are almost completely separated at the junction when the cyclotron diameter becomes shorter than 
the wire width (21 s (6) < N w ). 



B. The influence of disorder 



We now consider briefly the effect of disorder on the spin polarization (Fig|7J). An ensem- 
ble average is performed over 10 4 samples. The suppression of the polarization by disorder 
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becomes more prominent as the spin-orbit coupling becomes stronger. The mean free path 
of a 2DES in the tight-binding model is described by j^] 



3/2 ^ + 4^0 
48aV W 2 ■ 



(18) 



One can use this estimate to distinguish the ballistic regime from the diffusive one. For the 
present system, we obtain W ~ 1.53Vq for L m = 50a (indicated by an arrow in Fig. [7J). 
As seen in the figure, the sample size must be smaller than the mean free path in order to 
obtain high spin polarization. 
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FIG. 7: Conductance G and spin polarization \P\ as functions of the strength of disorder W for 
several spin-orbit couplings at E = — 3.8Vq. Ensemble average has been taken over 10 4 samples. For 
stronger spin-orbit coupling, the polarization becomes more sensitive to disorder. Arrow: crossover 
between ballistic and diffusive regimes [cf. Ea.l|18|)]. 



C. The influence of other channels 

Finally let us consider the system whose size is N w = 50a and Ni = 50a. Figure |H1 shows 
the dependence of the conductance and the spin polarization on the Fermi energy E for 
9 = 0.01-7T. In this energy region where the number of channels increases to values ranging 
from 5 to 10, several channels contribute to transport. While the spin polarization is reduced 
by channel mixing, it still stays higher than 10%. 
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FIG. 8: Conductance G and spin polarization |P| as functions of the Fermi energy E in the region 
of multi-channel transport. The system size is set to be N w = 50a and Ni = 50a, and the strength 
of spin-orbit coupling = 0.017T. The number of channel increases to values ranging from 5 to 10 
in this energy region. The polarization stays higher than 10%. 

V. SUMMARY 

In this paper we have investigated numerically spin polarized linear transport in a T- 
shaped conductor with Rashba spin-orbit coupling. We have considered single-channel 
transport and found that the spin polarization becomes almost perfect with a conductance 
close to Go for stronger spin-orbit coupling. Since the spin-orbit coupling can be regarded 
as the spin-dependent effective magnetic field the propagating electrons with certain 
initial spin states are deflected into one of the two output terminals by Lorentz force while 
those with opposite spin states are deflected into the other. The ratio between the cyclotron 
diameter in an effective field and the wire width depends on that between the 7r-phase spin 
precession length and the wire width. If the precession length becomes shorter than the 
wire width, the cyclotron diameter becomes shorter than the wire width so that electrons 
with antiparallel spins are almost completely separated at the junction and nearly 100% 
polarization can be achieved. 



With respect to the intrinsic spin Hall effect |27J, we should note the following point. 
The spin-orbit induced effective magnetic field can cause the anomalous Hall effect in fer- 
romagnets, where the spin direcion of conduction electrons is maintained by the interaction 
between localized spin states. However, it can not simply separate injected electrons with 
up and down spins parallel to the z-axis in the investigated system since the spin states of 




11 



propagating electrons are always changing due to the spin precession |25| . 

We have also investigated the effect of disorder on the polarization. Since the polarization 
becomes more sensitive to the disorder when the spin-orbit interaction is stronger, one needs 
to prepare the clean samples so that the system belongs to the ballistic regime to realize the 
perfect polarization of the current. 

In order to obtain information whether or not the predicted effects are observable in 
experiment, let us consider the parameters required for one of the favorable materials, InAs, 
with an effective mass m* = 0.039m (m is the free electron mass) and Rashba coupling 
a = 23.8 x 10 -12 eVm Let us assume that the width of the conduction band is A =1 eV. 
This gives for Vo = A/2Z = 125 meV {Z = 4 for square lattice), and for the tight-binding 
lattice parameter a = h/ '(2m*Vo) 1//2 ~ 2.8 nm. Using the above numerical value of a one 
obtains 9 = a/2Voa ~ 0.017T. This would reduce the polarization to about 10%, still a 
reasonable value for being observable in experiment. The crucial point, however, is the 
condition that the transport has to be in the single- channel regime. The wire width should 
be about 20 nm for single-channel transport when Fermi energy is of the order of 10 meV. 
In principle it is possible to fabricate such a narrow wire, but it makes the mean free path 
shorter and the effect of disorder may become critical. On the other hand, the wire width 
becomes 140 nm for the system with N w = 50a fSec lIV Cj) . This width of quantum wire 
can be easily fabricated than that for single- channel transport, which indicates that we can 
observe the spin polarization experimentally. 
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APPENDIX A: RECURSIVE GREEN FUNCTION METHOD FOR MULTI- 
TERMINAL GEOMETRY 



For numerical calculations, we apply the recursive Green function method 22j to the case 
of multi-terminal geometry. Let us consider three terminal geometry described in FigEl The 
central sample region (C) is attached to three probes (D, R and L). Each probe consists of 
the infinite ideal region and the sample one whose size is N p x (I = D,R and L). 



R 



D 



FIG. 9: Schematic draw of three terminal geometry. The shaded area represents a sample region 
where spin-orbit interaction is present. This sample region is connected to three electron reservoirs 
by ideal probes. 



The full Hamiltonian can be written down as 
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Ti. c is the Hamiltonian for the central sample region and can be described by 



n 



c 



c 



-V c 

v 21 



-vg 

ng 

-vg 



o 

-vg 



o 






\ 



(A2) 



\ 

where Hg denotes the Hamiltonian for the i-th slice along the y-direction and Vg the 
hopping term between the slice i and j. The hopping is restricted to nearest neighbours. 
V[ c is the hopping term between the central sample region and its neighbouring slices of 
three probes. 7"^ is the Hamiltonian for the set of i-th slices of three probes and can be 
described by 

( Hf ^ 
Hf 
Hf 

where Til (I = D,R and L) is the Hamiltonian for the i-th slice of the probe / (see FigfTUj) . 
Vij is the hopping term between the set of i-th slices and j-th ones. 



\ 



(A3) 



Np+l 1 i 1 1 



FIG. 10: Graphical interpretation for the Hamiltonian Hi. The i-th slices of three probes are 
gathered. 



Let us define some variables by following Ref . |22| as 

\ 



A 1 



(A4) 
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and 



uI = {<■■■ <l) » ( A5 ) 

where Af , . . . , X ! Nl are the eigenvalues of transfer matrix for the ideal region on the probe / 
and u[, . . . , wLj the eigenf unctions corresponding to A{, . . . , A^/ . We note that the incoming 
solutions and outgoing ones have the same form when there is no magnetic field. Then the 
effective Hamiltonian on site N p + 1 for the probe I can be written by 

VqF 1 , 



^Np+l 



(A6) 



with 



F 1 = U I A I (U I )~ 1 , (A7) 

where Vq is the hopping term in the ideal probe. Hn p +i can be obtained by exchanging H\ 
for H\ in Eq.flXHJ and setting i = N p + 1. 
We define corresponding Green function as 

Gr 



G 



1 



e-h 



'N P +1,N P +1 ^N P +1,N P 



GNp+1,1 Gn p +1,C ^ 



G 



N P ,N P +1 



Gr 



G 1 



Gr 



Gr 



G 



C,N V 



Gi,i 
Gc,i 



Gi,c 
Gc,c ) 



(A8) 



Gl,N p +l 
\ G C ,N P +1 

In principle, we can obtain Green function by the direct inversion of the matrix E — H. 
However, the number of numerical operations required for the inversion of the matrix in- 
creases as iV 3 with N the size of the matrix. We therefore need to reduce the size of matrix 
to inverse as small as possible. 

Since we put the effective Hamiltonian Hn p +i on the upper-left corner of the full Hamil- 
tonian fEq. (|Al|) ). the sub-matrix Gn p +i,n p +i contains all information of the transport. This 
sub-matrix can be calculated as 

1 



G 



Np+1,N P +1 



E — TLn„+i 



r 



Np+l 



where 



with 



i+1 
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V) > 
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The inverse matrix of E — Tt c can be calculated recursively |22| and each row of the ma- 
trix V[ c contains only one non-zero element corresponding to the matrix element between 
neighbouring slices, which greatly simplifies the calculation of V\. In the calculation of r i+1 
for i = 1,2, ...,N p and GjVp+i^+i) we need to inverse directly the matrix whose size is 
(N£ + N* + N%) x (N£ + N* + N%). This is much smaller than the original size olE-H. 
G n v +i,n v +i consists of 9 parts. 



/ qdd QDR qDL \ 



G 



N p +l,N p +l 



qRD qRR qRL 
qLD qLR qLL 



where G IJ (I, J = D, R and L) is the Green function from t 



matrix can be calculated easily from the Green function 



(A12) 

re probe J to /. The scattering 
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